Effective theory of magnetization plateaux in the Shastry-Sutherland lattice 
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We use the non-perturbative Contractor-Renormalization method (CORE) in order to derive an 
effective model for triplet excitations on the Shastry-Sutherland lattice. For strong enough magnetic 
fields, various magnetization plateaux are observed, e.g. at 1/8, 1/4, 1/3 of the saturation, as found 
experimentally in a related compound. Moreover, other stable plateaux are found at 1/9, 1/6 or 
2/9. We give a critical review of previous works and try to resolve some apparent inconsistencies 
between various theoretical approaches. 

PACS numbers: 75.60.Ej, 75.10.Jm 

The Shastry-Sutherland lattice []| and its realiza- 
tion in the material SrCu2 (603)2 have been attract- 
ing a lot of attention due to its fascinating behaviour 
in a magnetic field [2, y, Ij] , namely that magnetization 
plateaux have experimentally been observed for values of 
m = 1/8, 1/4, and 1/3 of the saturation value. The 
Shastry-Sutherland lattice, sketched in Fig. [TJa), is a 
two-dimensional (2d) Heiscnberg antifcrromagnctic spin- 
1/2 couplcd-dimcr system which Hamiltonian reads: 



H = J 



J2 s.-s, + j' J2 s..-s,-hj2s^. (1) 



<ij> 



<Si,i» 



Experiments with SrCu2(B03)2 indicate that the ratio 
between inter- and intra-dimer coupling should be close 
to J' /J sa 0.65 where the ground-state is exactly given by 
the product of singlets on J bonds [5| . In the presence of 
a finite magnetic field h, polarized triplets are created on 
the dimcr bonds so that low-energy properties can be de- 
scribed with an effective model of these "particles" : since 
these triplets arc hard-core bosons moving on an effective 
square lattice, they can typically exhibit compressible su- 
perfluid or incompressible Mott phases depending on the 
filhng [a, Q- In the original language, these two phases 
correspond respectively to absence or presence of magne- 
tization plateaux. 

In order to provide a simple picture, let us recall that 
the triplet hopping is strongly reduced on frustrated lat- 
tices, and therefore the physics is governed by effec- 
tive Coulomb repulsion, resulting in various insulating 
phases known as Wigner crystals. Following these pi- 
oneering works, other theoretical approaches have con- 
firmed the occurrence of several plateaux [5|, |7|. While 
all approaches agree to describe m = 1/2 or 1/3 plateaux, 
the situation is less clear at lower magnetization. 

Experimentally, because of accessible fields, the first 
plateaux were discovered at 1/8 and 1/4 Q], but recently, 
translation symmetry breaking has been observed above 
1/8 [SJ, |9[ as well as evidence of 1/6 platea u lldl . By us- 
ing torque measurements, Sebastian et al. [ll| have sug- 
gested additional magnetization plateaux for more "ex- 
otic" values hke 1/9, 1/7, 1/6, 1/5 and 2/9. Ahhough 



an agreement on these values has not been reached yet, 
these results are quite exciting and ask for a thorough 
theoretical analysis. 

Because of the lack of powerful numerical techniques 
or unbiased analytical tools to tackle 2d frustrated sys- 
tems, a promising approach consists in deriving an effec- 
tive hard-core bosonic model. Because one is interested 
in low-magnetization (i.e. low filling), long-range effec- 
tive interactions are crucially needed and can only be 
captured thanks to efficient algorithms. Recently, per- 
turbativc continuous unitary transformations up to high 
orders have provided such an effective model |12l | . In this 
letter, wc use the non-perturbative Contractor Rcnormal- 
ization (CORE) technique 13|, [ij] in order to derive an 
effective model for the polarized triplets. Then, in order 
to provide an unbiased analysis of this model, we solve it 
exactly on various clusters and by finite-size scaling anal- 
ysis, we predict the existence or not of some plateaux in 
the thermodynamical limit. 

Effective model - Similarly to perturbation theory [5|, 
|6[ , we base our approach on keeping only the singlet and 
polarized triplet states on each dimer [l5|. A crucial as- 
pect of the CORE technique is that it gives a cluster ex- 
pansion of the effective Hamiltonian Tieff- Basically, the 
amplitude of local processes can already be captured by 
solving a small finite cluster. As a consequence, the only 
approximation consists in truncating beyond a certain 
range of interactions [l3[ . In our study, we keep all pro- 
cesses that can appear on a 3 x 3 cluster (corresponding to 
18 original sites). Note that some elementary processes, 
as nearest-neighbor repulsion or chemical potential, are 
already well captured with a smaller 2x2 cluster, so that 
different truncations may give similar amplitudes. Nev- 
ertheless, long-range interactions are crucially needed to 
describe low-filling phases. Typically, Tieff contains of 
the order of lO'' terms, similar to what is obtained with 
high-order perturbation [I4I . 

Our CORE calculation gives access to the chemical po- 
tential ^ which corresponds to the spin gap. As sketched 
in Fig. [2lja), CORE results are in good agreement with 
4th order perturbation theory [5| up to ,/'/ J = 0.7. Be- 
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FIG. 1: (Color online) Shastry-Sutherland lattice and defi- 
nitions of some effective interactions: (a) 2-body interaction 
between a particle on the red dimer and one on the labelled 
dimers. (b) - (d) Typical correlated hopping processes. 



cause of the important role played by diagonal 2-body 
interactions, we compare them with high-order pertur- 
bation theory [l^l m Fig. [2Ia). Clearly, these interac- 
tions decrease with distance and are non-isotropic: for 
instance, V3 and V3' strongly differ. Generally, one also 
notices very good agreement between perturbation the- 
ory (third order [5|, IGJ or higher [I4I) and CORE results 
at least up to values of J' / J ~ 0.5. Beyond this value, 
other processes become dominant so that comparison be- 
comes more difficult to perform and we will restrict to 
J' / J < 0.5 in our study. 

Besides the discussed 2-body interaction terms, V^ ^, 
that corresponds to an attractive interaction of three 
aligned particles, has a sizeable contribution (-0.7) close 
to J' /J ~ 0.5, quite different from high-order pertur- 
bation by Dorier et al. (-0.4). Although being high- 
order processes in perturbation, we believe that inter- 
action terms involving three and more particles play an 
important role in the formation of plateaux. For exam- 
ple, in the one-dimensional realisation of the Shastry- 
Sutherland lattice, namely the orthogonal-dimer chain, 
diagonal interaction terms between several particles are 
responsible for the series of infinite plateaux [16| . 

The effective model also includes off-diagonal processes 
which account for hopping terms of bosons (i.e. triplets) 
from one dimer to another. In agreement with perturba- 
tion theory [6| , we find that simple one-particle hopping 
terms are negligible. On the contrary, correlated hopping 
(i.e. when one particle hops from one dimer to another 
in the presence of neighboring particle) is a dominant 
process and we illustrate the three most important such 
terms in Fig. [T] and give the corresponding amplitudes in 
Fig. [2jb). Even though they may appear small, corre- 
lated hopping terms are very important for the physics 
of a system as they favor supersolid phases \u\ . We note 
again a good agreement with perturbation theory up to 
J' / J w 0.5. Beyond that, as for diagonal terms, we ob- 
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FIG. 2: (Color online) Effective amplitudes, defined in Fig.[l] 
calculated with CORE (solid lines) and perturbation theory 
(dashed lines, see text): (a) ^ and 2-body diagonal terms; 
(b) correlated hoppings and high-order perturbation result 
for ta\h and taz from Dorier et al. (private communication). 



serve strong variations of the amplitudes that indicate 
the limit of validity of our CORE truncation. 

Validity of the CORE approach - It is known that 
the physics of the Shastry-Sutherland model changes 
from a dimer state to a 2d-like Heisenberg phase above 
(j'/j) ~ 0.70, possibly with an intermediate plaque- 
tte phase Q . Naturally, the question of the validity of 
our CORE approach emerges. In this regime, we observe 
that (i) basic processes amplitudes do not converge with 
different cluster sizes; (ii) ground-state has zero overlap 
with our subspace. Therefore, CORE procedure is only 
applicable below this critical value. 

Another useful tool to ascertain the validity of CORE 
approach is to compute exactly on small clusters the re- 
duced density matrix weights of retained states [1J| . A 
numerical analysis done on 16-site cluster confirms that 
the total weight of the 2 kept states exceeds 85% as long 
as J' / J < 0.65. This gives us confidence that effective 
interactions should decay fast enough so that our CORE 
procedure is accurate in this region. Because of the re- 
duced accuracy close to this J' /J value (which describes 
the SrCu2(B03)2 compound), we will restrict most of our 
findings to J' / J = 0.5 where various CORE truncations 
give similar models. Still, this value is reasonably close 
to the experimental one and no qualitative changes are 



expected in this region. In particular, the same magne- 
tization plateaux should occur. 

In all these effective models approaches, one must dis- 
tinguish the two steps: first, an effective Hamiltonian is 
derived; then, because it is still an interacting quantum 
problem, one needs to resort to an efficient technique to 
study it. Although it is a bosonic model, the presence 
of positive off-diagonal terms prohibits quantum Monte- 
Carlo. Possible alternatives are mean-field analysis [12 1 
which main drawback is to overestimate the tendency to 
form plateaux and also needs a careful numerical analysis 
on finite clusters. Given the various small amplitudes, we 
prefer not to make any further assumption and we pro- 
vide exact diagonalizations (ED) of these effective models 
on finite lattices. 

Simulations of the effective model - By solving ex- 
actly the effective models for various bosonic fillings, it is 
straightforward to construct the magnetization curve by 
a Legendre transform. A typical plot is given in Fig. [Sja) 
where we compare data obtained with the microscopic 
and effective models on A^ = 32 lattice. Because a given 
finite lattice cannot accommodate all magnetization val- 
ues (only multiples of 2/A^ are allowed), the magnetiza- 
tion curve presents many steps. Thus, we cannot con- 
clude yet about the existence of other plateaux, such as 
1/6, and one needs to do a careful finite-size extrapola- 
tion to get information on the thermodynamical limit. 

Nevertheless, for this system size, Fig.[3l^b) shows that 
our effective Hamiltonian is extremely accurate, as it re- 
markably coincides with exact values (both for plateaux 
widths and locations) at least for magnetization m < 0.5. 
Moreover, this accuracy is excellent up to J' /J ~ 0.4 and 
we only observe a small disagreement for m — 1/2 and 
J' /J = 0.5. A huge advantage of our effective model 
is that, although it involves many terms (basically, all 
terms up to 9-body), due to Hilbert space reduction, we 
are able to solve systems twice as large as for the micro- 
scopic model, which makes feasible a finite-size scaling. 
We have therefore solved Tios on a. N = 64 lattice and a 
its magnetization curve is given in Fig. [3Ka). Naturally, 
there are twice as many finite-size steps but the main 
message is that to = 1/4 and to = 1/8 plateaux do not 
change, both in sizes and locations. This strongly sug- 
gests that our model does exhibit magnetization plateaux 
in the thermodynamical limit. Because we have chosen 
particular clusters, there is also the possibility to have 
more stable plateaux close to these fractions, and for in- 
stance, some plateaux that we find could be unstable 
towards phase separation. 

Other possible fractions can be investigated by per- 
forming similar calculations with other square clusters, 
like A^ = 36 and 72, that can also accommodate to = 1/9, 
1/6, 2/9 and 1/3. Data are shown in Fig.jSjc-d) and con- 
firm (i) the accuracy of our CORE effective model and 
(ii) the stability of some plateaux when the system size is 
doubled, strongly indicating that they do persist in the 
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FIG. 3: (Color online) (a) Magnetization curve obtained 
with microscopic and CORE models on A = 32 lattice and 
CORE results on 64-site for J' / J — 0.5. (b) Phase diagram 
for A = 32 as a function of J' /J and magnetic field h/J: 
CORE results (lines) agree with ED (symbols) for locations 
of m = 1/2, 1/4 and 1/8 plateaux which are allowed on this 
cluster, (c)-(d) Same as (a)-(b) for A = 36 and 72. On these 
clusters, m = 1/2, 1/3, 1/6, 1/9 or 2/9 plateaux are allowed. 



thermodynamical limit. 

In order to give a general view of various fractions, 
we now turn to a more systematic study of our effective 
Hamiltonian, restricted to the typical J' /J = 0.5 value, 
on various square or rectangular clusters. By computing 
the magnetization curves on several square lattices, we 
can perform a finite-size scaling of the plateaux widths 
(see Fig. 21). Note that we restrict our calculations to clus- 
ters that can accommodate a given insulating phase, i.e. 
arc not frustrated, according to the known patterns [3|,|5[. 
Since we also have access to density-density correlations, 
we can also confirm these patterns (data not shown): for 
instance, at m = 1/3 (resp. m = 1/4), the plateau is 
formed by filled diagonal stripes separated by two (resp. 
three) empty ones. 

Our scaling provides clear evidence that there are large 
plateaux for to = 1/3 and 1/4 (besides m = and 1/2) 
in the thermodynamical limit [l8|. Moreover, within the 
very good accuracy of our effective model, finite-size scal- 
ing indicates smaller, but stable, plateaux for 2/9, 1/6, 
1/8 and 1/9 of typical sizes between 0.01 and 0.05 J. 

Discussion - By allowing for inhomogeneous pat- 
terns, a variety of fractions have been found thcorcti- 
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FIG. 4: (Color online) Finite-size scaling of the plateaux 
widths for various m < 1/2 from microscopic and CORE 
models (respectively open and filled symbols) for J' /J = 0.5 
on clusters with 24, 32, 36, 48, 64, 72 and 96 sites (square 
and rectangular). Note the excellent agreement between both 
models solved on same lattices. 



cally ll|, [l^. Experimentally, there is recent evidence 
that translation symmetry is broken for several other 
magnetization fractions 1Q|, |ll|. Clearly, given the small 
energy scales that stabilize one fraction with respect to 
another, the accuracy of the calculation is crucial. In 
that sense, systematic perturbative [l^l or CORE expan- 
sion (as in this Letter) are promising since they are con- 
trolled techniques. For instance, our relative error on the 
ground-state energy for TV = 32 and m = 1/8 is 0.3%. 

However, from these effective models, a second step 
is to solve them in the most unbiased way: here we 
have combined exact diagonalization and finite-size scal- 
ing and we find highly stable plateaux for 1/2, 1/3, 1/4, 
but also smaller plateaux for 2/9, 1/6, 1/8 and 1/9. 
While some of these fractions coincide with [l^ , we also 
have some differences that ask for a clarification. 

A first discrepancy is the absence in our data of 
plateaux for m = 2/15, observed in |12| or possible 
other fractions such as 1/7 or 1/5 found in [ll|. Clearly, 
because our exact simulations are restricted to large but 
finite clusters, we cannot perform finite-size scaling for 
some fractions (as 2/15 that would require large unit 
cells). By choosing adequate shapes or boundary con- 
ditions, it might be possible to investigate some of these 
other fractions. 

A more crucial issue deals with m = 1/4 and 1/8 
plateaux that are stable in our calculations, and also 
found experimentally [3|, but absent in [13|. A possi- 
ble explanation could be that these fractions are unsta- 
ble towards phase separation but, on our finite clusters, 
this phenomenon cannot be observed. However, these 
plateaux widths are almost constant when the system 
size is doubled (see Fig. 2]) so that we do strongly believe 



that they persist in the thermodynamical limit. At this 
point, let us recall that our effective model is different 
from [l^l, in particular for interactions beyond 2-body: 
for instance, if we only consider 2-body diagonal terms, 
then the m = 1/4 plateau becomes unstable, in agree- 
ment with [l^]. Therefore, we think that many-body 
terms ask for a careful analysis. 

As a conclusion, we are convinced that reliable effec- 
tive hamiltonians are crucial to understand the very rich 
low-magnetization properties of the Shastry-Sutherland 
system. With the CORE technique, we have derived such 
an effective model and, with finite-size scaling analysis, 
we provide a microscopic origin of the experimentally ob- 
served 1/3, 1/4 and 1/8 plateaux, but we also confirm 
the possibility of other plateaux at 1/9, 1/6 and 2/9, as 
found in recent related studies 10|, [llj, [l^ . The remain- 
ing discrepancies between these approaches call for a sys- 
tematic unbiased study, that would combine mean-field 
ideas with exact diagonalizations. 
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